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1.  Introduction 


A  number  of  authors  hava  reviewed  the  structure  of  the  wake 
generated  by  relative  notion  between  an  object  and  a  collisionless 
plasa».t*"®3  These  reviews  have  conaidered  the  wake  structure  extending 
from  the  near  wake  to  the  far  wake  for  various  ranges  of  body  size  R0 
relative  to  the  ambient  plasma  Debye  length  Xq  and  for  various 
conditions  of  electrical  potential  on  the  body  relative  to  the  ambient 
plasma. 

This  study  focuses  on  the  structure  of  the  ion  front  in  the  near 
wake  behind  a  semi-infinite  plate,  and  in  particular  on  the  effect  of 
plate  potential  on  the  structure  of  the  ion  front.  The  present 
analysis,  which  is  relevant  to  the  wake  structure  of  a  large  body 
traversing  the  ionosphere,  is  motivated  by  the  experiments  on  wakes 
performed  in  a  large  vacuum  tank^.  Fig.  1  depicts  the  experimental 
situation.  The  measured  results  obtained  by  means  of  a  differential  ion 
flux  probe  consist  of  ion  density  profiles  as  a  function  of  x  at  a  given 
downstream  position.  The  plate  is  biased  to  a  potential  ^p  relative  to 
the  tank  wall,  and  ^p  may  differ  from  the  potential  of  the  plasma  in 
the  neighborhood  of  the  plate.  The  effect  of  lp  (or  lp-^a)  on  the 
aiot ion  and  structure  of  the  ion  front  in  the  wake  behind  the  plate  is 
determ i ned . 

The  analysis,  exploiting  the  approximate  relationship  between  the 
problems  of  a  mesosonic  wake  and  the  expansion  of  plasma  into  vacuum 


considers  first  the  effect  of  an  applied  electric  field  on  the  initial 
piassia  expansion.  Specific  calculations  of  wake  structure,  including 
the  effect  of  an  applied  potential  and  of  ion  temperature,  are  given  and 
compared  with  experiment  in  Section  3. 


2.  Effect  of  Potential  on  Expansion  into  Vacuum 


For  zero  bias,  the  electric  field  at  the  ion  front  xp(z)  is  given 
approximately  byt7^ 


a!  ,  2 

M  1  ♦  awt 


1  “  (1_  ^  7 )  rrsit] 


(1) 


t  .  z/M 


with  the  adjustable  parameter  a,  e  »  2.718,  and  ions  travel  with 

velocity  VQ  in  the  positive  Z  direction.  Here  q  and  M  are  the  ion 
charge  and  mass,  respectively,  0  the  electron  temperature  and 

2^/2 

Nq 

«  *  jy  the  ion  plasma  frequency  for  ambient  density  N.  The  semi- 
.  o 

empirical  formula  (1)  was  constructed  to  reproduce,  independent  of  a, 
certain  theoretically  known  properties  of  expansion  of  a  plasma  half 
space  into  vacuum,  namely  the  asymptotic  dependence  of  the  ion  front  on 
t  for  wt  >>  it8!,  and  the  value  of  E  at  t  =  0,  namely  E0  = 
(2H6/eQe)  The  position  end  ve  locity  of  the  front  calculated  from 

eq.  (1)  are  in  good  agreement  with  the  measurements  of  Wright,  Stone  and 


Let  us  now  consider  how  the  expansion  into  vacuum  is  modified  by  an 
applied  field  Ej  in  the  vacuum  region.  Initially  the  ion  density  is  N 
(uniform)  for  x  <  0,  and  zero  for  x  >0.  The  potential  is  determined 
from 


V2*  a  Nq  (l-exp  (|^)  )  X  <  0  j 
a  Nq  exp(^)  x  >  0 


(2) 


and  satisfies  the  conditions  ^(x=-“)  =  0,  ^(x=«)  =  -* ,  E(x=-«)  =  0, 
E(xa»)  =  E. .  The  first  integral  of  (2)  over  the  plasma  half  space  from 

x  =  -•  to  x=0  yields 

'f  E02  =  NS[.xp(q*o/0)  -  1  -  q  ^/«]  (3) 


where  a  subscript  o  denotes  the  value  at  the  plasma  vacuum  interface. 
Integration  over  the  vacuum  region  yields 

1?  (Eo2  -  Ej2)  a  N*  exp(q 4J6) 

For  Ej  *  0,  *e  a  ^Qo  and  Eq  a  Eq<>  where 


(4) 


00 


a  -1 


-f  E  2  =  Nfle"1 
2  oo 


(5) 


are  the  results  obtained  by  Crow,  Auer  and  Allen^. 
For  E.  #  0, 

i  ' 


(6) 


The  applied  field  is  the  dominant  component  of  the  initial  field  if 

A 

e0Ej  /2N0  >  1.  The  applied  field  is,  of  course,  a  nebulous  quantity;  in 
particular,  it  is  not  unambiguously  related  to  the  value  of  the  bias 
potential  on  an  object  over  which  a  plasma  is  flowing.  In  order  of 
magnitude  however,  E;  should  be  comparable  to  that  of  the  sheath 
electric  field  which  deflects  the  streaming  ions  before  they  enter  the 
near  wake  downstream  from  the  object;  E;  ~  l0/R  where  is  the 

potential  on  the  object  relative  to  that  in  the  remote  upstream  plasma 
and  R  is  a  length  scale  for  potential  variation  within  the  sheath. 
Rather  than  attempting  to  further  pursue  the  analogy  between  the  plasma 
expansion  and  the  wake  problems,  we  turn  now  to  consideration  of  the  two 
dimensional  steady  wake. 

3 .  Qualitative  Effect  of  Potential  on  a  Two-Dimensional  Steady  State 

Wake 

The  expansion  of  a  semi  - i nf i ni te  plasma  into  vacuum  corresponds  only 
approximately  to  the  wake  problem.  In  particular,  the  free  expansion 
analogue  does  not  take  into  account  the  fixed  potential  on  the  wake- 
forming  plate,  nor  does  it  clearly  suggest  how  to  account  for  variations 


IV, 


in  plate  potential.  A  complete  theory  would  be  based  on  the  steady 
state  solutions  of  the  Poisson  and  ion  trajectory  equations  in  at  least 
two  dimensions. 

A  primary  effect  neglected  in  the  free  expansion  model  is  the 
transverse  acceleration  experienced  by  ions  that  pass  from  upstream 
(z  <  0)  through  the  sheath  surrounding  the  tip  of  the  plate  before 
entering  the  wake  region  (z  >  0,  x  >  0)  (see  Fig.  1).  The  sheath  should 
extend  in  the  upstream  direction  no  more  than  a  few  Xq,  for  £  100. 
In  these  circumstances  the  approximate  effect  of  the  sheath  on  fast  ions 
(1/2  MV0^  »  Ip)  can  be  described  in  terms  of  an  impulsive  change  in 
transverse  velocity,  Vx: 


M  AVX  2  qEAt 


(9) 


^D1  R_ 
R  V 


where  R  is  the  range  of  the  force,  i.e. 


AV  q* 

_*  ~  rP 

Vo  'MV2 

O 


(10) 


For  the  following  parameters,  corresponding  to  conditions  in  Wright’s 


exper iments 


[6] 


4  =  -3  volts 

rP 


|  MVq2  =  18  volts  (N2*) 


5 


At  a  point  z  =  22.4  cm  downstream,  Ax  (est)  =  y-  AV^  -  2  cm. 

o 

Assuming  that  the  effect  of  the  finite  plate  potential  on  the  wake 
is  to  deliver  impulse  to  the  ions,  this  means  that  the  position  xp  of 
the  ion  front  will  have  shifted  by  an  amount  Ax  relative  to  that  for  an 
unbiased  plate.  The  shift  in  experimental  values  of  xp  between  ^  =  0 
and  Ip  =  -3  volts  is 

Ax  (exp)  =  .65  cm. 

Although  of  the  right  order  of  magnitude,  the  estimated  Ax  (est.)  is  3 
times  Ax  (obs)  even  before  taking  account  of  the  effect  of  potential  in 
the  downstream  direction.  In  any  case,  the  impulse  approximation  given 
here  leads  to  a  displacement  linear  in  applied  potential,  a  conclusion 
not  inconsistent  with  experimental  results,  at  least  for  negative 
potentials. 

For  positive  potentials,  the  effect  of  potential  appears  to  be 
smaller  than  for  negative  potentials,  except  for  an  abrupt  change 
between  ^  =  0  and  ^  *  +1  volt.  A  smaller  effect  could  be  expected  on 
the  grounds  that  in  the  neighborhood  of  the  plate  the  plasma  will  be 
only  a  few  $  below  ^  for  ^  >  0. 
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3.  Two-Dimensional  Calculations  of  Near  Wake  Structure 


The  experiment  was  simulated  numerically  in  two  dimensions  using  a 
finite  element  computer  code.  The  code  solves  the  Poisson-Vlasov 
equations  in  steady  state.  Potentials,  ion  trajectories,  and  ion 
densities  were  calculated  self-consi stent  I y  on  the  finite  element  mesh 
shown  in  Figure  2.  The  boundary  conditions  were  zero  potential  at  the 
top  (ambient  plasma)  and  left  (input)  boundaries,  specified  (imposed 
bias)  potential  on  the  metal  plate,  zero  normal  component  of  electric 
field  on  the  right  (exit)  boundary,  and  =2  exp  (f/O)  on  the  lower 
(vacuum)  boundary.  For  the  finite  element  potential  solution,  the 
potential  was  assumed  to  vary  bi linearly  in  each  element.  For 
calculating  trajectories,  however,  a  quadratic  form  was  used  to 
interpolate  in  the  direction  approximately  normal  to  the  ion  front. 
This  was  necessary  because  the  discontinuous  electric  fields  obtained 
from  linear  interpolation  caused  bunching  of  ion  trajectories,  resulting 
in  striations  in  the  downstream  ion  density. 

The  ion  density  was  calculated  using  a  variant  of  the  "waterbag" 
method [10,11]  This  method  takes  advantage  of  the  fact  that,  if 
trajectories  do  not  cross,  the  current  contained  between  two  trajectory 
paths  (in  two  dimensions,  or  in  a  "flux  tube"  in  three  dimensions)  is 
constant.  Denoting  by  u  the  dimension  along  the  trajectories  and  by  w 
the  dimension  normal  to  the  trajectories,  we  have 


IAt  =  nvAwAt  =  nAwAu 


where  n  and  v  are  the  local  ion  density  and  velocity.  It  follows  that 
the  density  ratio  for  quadrilaterals  formed  by  adjacent  time  points  on 
adjacent  trajectories  (Figure  3)  is  the  inverse  ratio  of  their  areas. 
Calculations  were  performed  with  initial  (left  boundary  of  grid)  values 
of  0.1  and  0.2  Debye  lengths  for  Aw  and  Au  respectively,  with  constant 
timesteps  At.  For  the  finite  temperature  case,  the  ions  were  divided 
into  ten  subspecies,  each  representing  a  different  value  of  transverse 
velocity,  and  the  resultant  subspecies  densities  were  added  together. 
This  technique  for  calculating  ion  densities  resulted  in  very  little 
numerical  noise  and  no  code  diffusion,  and  proved  far  superior  to 
particle-in-cel  I  methods  for  cases  such  as  this  where  a  sharp  ion 
density  boundary  exists. 

The  electron  density  was  assumed  to  satisfy  the  "barometric  law" 

ne/N  =  exp(q^/0) 

for  negative  potentials,  and 

ne/N  =  1  ♦  q?/0 

for  positive  potentials,  where  N  and  9  are  the  ambient  plasma  density 
(10-5  cnT^)  and  electron  temperature  (0.27  eV) ,  and  f  is  the  local 


potential.  Ion  densities  and  potentials  were  iterated  until  neither 
changed  more  than  one  percent. 


Figure  4  shows  the  electrostatic  potentials  calculated  for  the  zero 
temperature,  zero  bias  case.  There  are  no  electric  fields  seen  upstream 
of  the  plate.  Downstream,  beyond  a  transition  region  of  about  one  Debye 
length,  the  equ i potent i a  I  lines  are  nearly  parallel  to  the  ion  flow 
direction,  and  equivalence  to  the  one-dimensional  time-dependent 
problem^  is  nearly  exact.  (Note  that  the  influence  of  the  plate 
extends  several  Debye  lengths  into  the  vacuum  region.)  Electrostatic 
potential  for  the  case  of  the  plate  biased  to  -106  is  shown  in  Figure  5. 
In  this  case,  the  ions  see  very  strong  downward  fields  as  they  pass  over 
the  plate,  providing  the  impulse  described  in  Equation  (10).  Moreover, 
the  influence  of  the  plate  potential  is  strongly  felt  by  the  ions  for 
about  ten  Debye  lengths  downstream,  providing  additional  downward  force 
on  the  ion  front. 

Ion  density  profiles  calculated  15  Debye  lengths  downstream  of  the 
plate  are  shown  in  Figure  6  for  cold  ions,  and  in  Figure  7  for  a 
transverse  ion  temperature,  9\  of  0.096  6  (300  degrees  Kelvin).  The  two 
curves  in  each  figure  are  for  plate  potentials  of  zero  and  -10  8.  The 
increase  in  front  motion,  (about  1.8  Debye  lengths),  corresponds  closely 
with  the  2  cm  obtained  from  Equation  (10).  The  finite  temperature 
profiles  are  monotonic  and  considerably  broader  than  the  cold  ion 
profiles.  The  angular  spread  of  the  ion  front  is  a  few  times 
(20j/MVo^),  as  expected.  The  only  non-monotonic  behavior  is  for  cold 
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ions  with  zero  potential  on  the  plate.  Figure  8(a-d)  show  the 
calculated  densities  throughout  the  region  for  the  four  cases. 

The  calculated  finite  temperature  ion  density  is  compared  with  the 
laboratory  measurement  in  Figure  9.  The  general  front  amplitude  and 
motion  agrees  well  with  experiment.  However,  the  calculation  has  a 
broader  front,  while  the  experiment  evidences  an  overshoot  and 
undershoot. 

4.  Discussion 

The  calculated  and  experimental  results  show  that  wake  closure  is 
well  described  by  the  acceleration  of  ions  in  the  plasma  steady-state 
electric  field.  However,  the  electrostatic  influence  of  the  object 
creating  the  wake  leads  to  significant  differences  between  wake  closure 
and  the  expansion  of  a  plasma  into  vacuum. 

For  zero  ion  temperature  and  zero  bias  on  the  plate,  the  calculated 
front  position  (0.50Ap)  is  less  than  theoretically  predicted  (0.72Aq). 
This  is  due  to  the  reduction  of  the  transverse  electric  field  at  the 
front  by  the  influence  of  the  conducting  plate  for  several  Xq 
downstream.  The  experimentally  measured  front  position  (1.5\[))  is 
appreciably  greater  than  predicted. 

The  application  of  -100  bias  to  the  plate  dramatically  changes  the 
electrostatic  potential  structure.  The  shift  in  front  position 
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predicted  by  the  crude  "impulse"  theory  (l.lXp)  is  in  qualitative 
agreement  with  calculation  (1.8Xq).  The  experimental  shift  (0.5Xq)  is 
smaller,  but  the  resultant  front  positions  (2Xp  experiment;  2.3Xq 
calculation)  agree  well.  Experiment  and  calculation  also  show  excellent 
agreement  as  to  the  magnitude  of  density  change  at  the  ion  front. 

The  calculations  produce  the  expected  thermal  spread  in  the  front. 
The  experimental  thermal  spread  is  far  less,  suggesting  that  the 
transverse  ion  temperature  is  welt  below  300K.  (Had  we  shown  the  zero 
temperature  calculation  on  Figure  9,  apparent  agreement  would  have  been 
near  perfect.)  However,  an  ion  temperature  of  300K  produces  no 
observable  change  (relative  to  zero  temperature  ions)  in  the 
electrostatic  potential.  We  believe  that  the  reason  for  lack  of  thermal 
spread  in  the  experimental  wake-front  is  that  the  ions  originate  from  a 
small  source,  so  that  they  are  already  "velocity-selected"  when  they 
enter  the  wake  region. 

Expansion  of  plasma  into  vacuum  is  a  good  approximate  model  for 
filling  in  a  plasma  wake.  However,  multi-dimensional  effects  due  to  the 
influence  of  the  body  creating  the  wake  on  the  electrostatic  potential 
structure  cause  noticeable  departures  from  this  model  when  the  body  is 
at  plasma  potential,  and  substantial  departures  when  the  body  departs 
from  plasma  potential. 
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FIGURES 


Figure  1.  Schematic  illustration  of  experiment  to  measure  ion  density 
in  wake  behind  a  biased  plate. 

Figure  2.  Finite  element  mesh  used  to  calculate  potentials  and  ion 
trajectories  in  the  wake  behind  a  biased  plate.  (Distances 
are  in  Debye  lengths.) 

Figure  3.  Schematic  of  method  for  determining  ion  densities. 

Quadrilaterals  are  drawn  between  adjacent  time  points  of 
adjacent  trajectories.  The  density  ratio  is  the  inverse  of 
the  area  ratio. 

Figure  4.  Electrostatic  potential  contours  for  the  zero  bias 

potential  case  at  zero  transverse  ion  temperature.  Contour 
levels  are  at  intervals  of  half  the  electron  temperature. 
The  same  plot  for  finite  ion  temperature  is 
indistinguishable  from  this  one. 

Figure  5.  Electrostatic  potentials  contours  for  bias  potential  of  - 
106  at  finite  (300K)  transverse  ion  temperature.  Contour 
levels  are  at  intervals  of  the  electron  temperature  The 
same  plot  for  zero  ion  temperature  is  indistinguishable 
from  this  one. 


Figure  6. 


Figure  7. 


Figure  8. 


Figure  0. 


Ion  density  profiles  for  zero  ion  temperature  calculated  15 
Debye  lengths  downstream  of  the  piste.  The  two  curves  are 
for  plate  bias  potentials  of  zero  (showing  a  peak  and  sharp 
drop  at  -0.5  Debye  lengths)  and  -108  (showing  a  sharp  drop 
at  -2.3  Debye  lengths). 

Ion  density  profiles  for  finite  (300K)  ion  temperature 
calculated  15  Debye  lengths  downstream  of  the  plate.  The 
two  curves  are  for  plate  bias  potentials  of  zero  (showing  a 
broadened  front  from  *0.1  to  -1.3  Debye  lengths)  and  -108 
(showing  a  broadened  front  from  -1  to  -3  Debye  lengths). 

Ion  densities  for 

(a)  zero  temperature,  zero  bias; 

(b)  finite  (300K)  temperature,  zero  bias; 

(c)  zero  temperature,  -108  bias, 

(d)  finite  (300K)  temperature,  -108  bias. 

The  spacing  between  contour  levels  is  approximately  0.06, 
with  several  contour  lines  on  the  discontinuity  ((a)  and 
(c))  or  in  the  shaded  region  ((b)  and  (d)). 

Measured  ion  density  profile  for  -3V  bias,  together  with 
the  calculated  curve  for  300K  ion  temperature  and  -108 
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